Near-infrared emission from sublimating dust in collisionally active debris disks

Hot exozodiacal dust is thought to be responsible for excess near-infrared (NIR) emission emanating from the innermost parts of some debris disks. The origin of this dust, however, is still a matter of debate. We test whether hot exozodiacal dust can be supplied from an exterior parent belt by Poynting-Robertson (P-R) drag, paying special attention to the pile-up of dust that occurs due to the interplay of P-R drag and dust sublimation. Specifically, we investigate whether pile-ups still occur when collisions are taken into account, and if they can explain the observed NIR excess. We compute the steady-state distribution of dust in the inner disk by solving the continuity equation. First, we derive an analytical solution under a number of simplifying assumptions. Second, we develop a numerical debris disk model that for the first time treats the complex interaction of collisions, P-R drag, and sublimation in a self-consistent way. From the resulting dust distributions we generate thermal emission spectra and compare these to observed excess NIR fluxes. We confirm that P-R drag always supplies a small amount of dust to the sublimation zone, but find that a fully consistent treatment yields a maximum amount of dust that is about 7 times lower than that given by analytical estimates. The NIR excess due this material is much smaller (<10^-3 for A-type stars with parent belts at >1 AU) than the values derived from interferometric observations (~10^-2). Pile-up of dust still occurs when collisions are considered, but its effect on the NIR flux is insignificant. Finally, the cross-section in the innermost regions is clearly dominated by barely bound grains.Comment: 18 pages, 10 figures, A&A accepte

Laboratory-based grain-shape models for simulating dust infrared spectra

Analysis of thermal dust emission spectra for dust mineralogy and physical grain properties depends on laboratory-measured or calculated comparison spectra. Often, the agreement between these two kinds of spectra is not satisfactory because of the strong influence of the grain morphology on the spectra. We investigate the ability of the statistical light-scattering model with a distribution of form factors (DFF model) to reproduce experimentally measured infrared extinction spectra for particles that are small compared to the wavelength. We take advantage of new experimental spectra measured for free particles dispersed in air with accompanying information on the grain morphology. For the calculations, we used DFFs that were derived for aggregates of spherical grains, as well as for compact grain shapes corresponding to Gaussian random spheres. Irregular particle shapes require a DFF similar to that of a Gaussian random sphere with sigma=0.3, whereas roundish grain shapes are best fitted with that of a fractal aggregate of a fractal dimension 2.4-1.8. In addition we used a fitting algorithm to obtain the best-fit DFFs for the various laboratory samples. In this way we can independently derive information on the shape of the grains from their infrared spectra. For anisotropic materials, different DFFs are needed for the different crystallographic axes. This is due to a theoretical problem, which is inherent to all models that are simply averaging the contributions of the crystallographic directions.Comment: 8 pages, 8 figures, accepted by Astronomy and Astrophysic

Dust-grain processing in circumbinary discs around evolved binaries. The RV Tauri spectral twins RU Cen and AC Her

Context: We study the structure and evolution of circumstellar discs around evolved binaries and their impact on the evolution of the central system. Aims: To study in detail the binary nature of RUCen and ACHer, as well as the structure and mineralogy of the circumstellar environment. Methods: We combine multi-wavelength observations with a 2D radiative transfer study. Our radial velocity program studies the central stars, while our Spitzer spectra and broad-band SEDs are used to constrain mineralogy, grain sizes and physical structure of the circumstellar environment. Results: We determine the orbital elements of RUCen showing that the orbit is highly eccentric with a rather long period of 1500 days. The infrared spectra of both objects are very similar and the spectral dust features are dominated by Mg-rich crystalline silicates. The small peak-to-continuum ratios are interpreted as being due to large grains. Our model contains two components with a cold midplain dominated by large grains, and the near- and mid-IR which is dominated by the emission of smaller silicates. The infrared excess is well modelled assuming a hydrostatic passive irradiated disc. The profile-fitting of the dust resonances shows that the grains must be very irregular. Conclusions: These two prototypical RVTauri pulsators with circumstellar dust are binaries where the dust is trapped in a stable disc. The mineralogy and grain sizes show that the dust is highly processed, both in crystallinity and grain size. The cool crystals show that either radial mixing is very efficient and/or that the thermal history at grain formation has been very different from that in outflows. The physical processes governing the structure of these discs are similar to those observed in protoplanetary discs around young stellar objects.Comment: 11 pages, 12 figures, accepted for publication by A&

The perturbed sublimation rim of the dust disk around the post-AGB binary IRAS08544-4431

Context: Post-Asymptotic Giant Branch (AGB) binaries are surrounded by stable dusty and gaseous disks similar to the ones around young stellar objects. Whereas significant effort is spent on modeling observations of disks around young stellar objects, the disks around post-AGB binaries receive significantly less attention, even though they pose significant constraints on theories of disk physics and binary evolution. Aims: We want to examine the structure of and phenomena at play in circumbinary disks around post-AGB stars. We continue the analysis of our near-infrared interferometric image of the inner rim of the circumbinary disk around IRAS08544-4431. We want to understand the physics governing this inner disk rim. Methods: We use a radiative transfer model of a dusty disk to reproduce simultaneously the photometry as well as the near-infrared interferometric dataset on IRAS08544-4431. The model assumes hydrostatic equilibrium and takes dust settling self-consistently into account. Results: The best-fit radiative transfer model shows excellent agreement with the spectral energy distribution up to mm wavelengths as well as with the PIONIER visibility data. It requires a rounded inner rim structure, starting at a radius of 8.25 au. However, the model does not fully reproduce the detected over-resolved flux nor the azimuthal flux distribution of the inner rim. While the asymmetric inner disk rim structure is likely to be the consequence of disk-binary interactions, the origin of the additional over-resolved flux remains unclear. Conclusions: As in young stellar objects, the disk inner rim of IRAS08544-4431 is ruled by dust sublimation physics. Additional observations are needed to understand the origin of the extended flux and the azimuthal perturbation at the inner rim of the disk.Comment: Accepted for publication in A&A, 13 figures, 13 page

The evolved circumbinary disk of AC Her: a radiative transfer, interferometric and mineralogical study

We aim to constrain the structure of the circumstellar material around the post-AGB binary and RV Tauri pulsator AC Her. We want to constrain the spatial distribution of the amorphous as well as of the crystalline dust. We present very high-quality mid-IR interferometric data that were obtained with MIDI/VLTI. We analyse the MIDI data and the full SED, using the MCMax radiative transfer code, to find a good structure model of AC Her's circumbinary disk. We include a grain size distribution and midplane settling of dust self-consistently. The spatial distribution of crystalline forsterite in the disk is investigated with the mid-IR features, the 69~$\mu$m band and the 11.3~$\mu$m signatures in the interferometric data. All the data are well fitted. The inclination and position angle of the disk are well determined at i=50+-8 and PA=305+-10. We firmly establish that the inner disk radius is about an order of magnitude larger than the dust sublimation radius. Significant grain growth has occurred, with mm-sized grains being settled to the midplane of the disk. A large dust mass is needed to fit the sub-mm fluxes. By assuming {\alpha}=0.01, a good fit is obtained with a small grain size power law index of 3.25, combined with a small gas/dust ratio <10. The resulting gas mass is compatible with recent estimates employing direct gas diagnostics. The spatial distribution of the forsterite is different from the amorphous dust, as more warm forsterite is needed in the surface layers of the inner disk. The disk in AC Her is very evolved, with its small gas/dust ratio and large inner hole. Mid-IR interferometry offers unique constraints, complementary to mid-IR features, for studying the mineralogy in disks. A better uv coverage is needed to constrain in detail the distribution of the crystalline forsterite in AC Her, but we find strong similarities with the protoplanetary disk HD100546.Comment: update with final version published in A&

Dusty tails of evaporating exoplanets. II. Physical modelling of the KIC 12557548b light curve

Evaporating rocky exoplanets, such as KIC 12557548b, eject large amounts of dust grains, which can trail the planet in a comet-like tail. When such objects occult their host star, the resulting transit signal contains information about the dust in the tail. We aim to use the detailed shape of the Kepler light curve of KIC 12557548b to constrain the size and composition of the dust grains that make up the tail, as well as the mass loss rate of the planet. Using a self-consistent numerical model of the dust dynamics and sublimation, we calculate the shape of the tail by following dust grains from their ejection from the planet to their destruction due to sublimation. From this dust cloud shape, we generate synthetic light curves (incorporating the effects of extinction and angle-dependent scattering), which are then compared with the phase-folded Kepler light curve. We explore the free-parameter space thoroughly using a Markov chain Monte Carlo method. Our physics-based model is capable of reproducing the observed light curve in detail. Good fits are found for initial grain sizes between 0.2 and 5.6 micron and dust mass loss rates of 0.6 to 15.6 M_earth/Gyr (2-sigma ranges). We find that only certain combinations of material parameters yield the correct tail length. These constraints are consistent with dust made of corundum (Al2O3), but do not agree with a range of carbonaceous, silicate, or iron compositions. Using a detailed, physically motivated model, it is possible to constrain the composition of the dust in the tails of evaporating rocky exoplanets. This provides a unique opportunity to probe to interior composition of the smallest known exoplanets.Comment: 18 pages, 11 figures, A&A accepte

Chromatic transit light curves of disintegrating rocky planets

Context. Kepler observations have revealed a class of short period exoplanets, of which Kepler-1520 b is the prototype, which have comet-like dust tails thought to be the result of small, rocky planets losing mass. The shape and chromaticity of the transits constrain the properties of the dust particles originating from the planet's surface, offering a unique opportunity to probe the composition and geophysics of rocky exoplanets. Aims. We aim to approximate the average Kepler long-cadence light curve of Kepler-1520 b and investigate how the optical thickness and transit cross-section of a general dust tail can affect the observed wavelength dependence and depth of transit light curves. Methods. We developed a new 3D model that ejects sublimating particles from the planet surface to build up a dust tail, assuming it to be optically thin, and used 3D radiative transfer computations that fully treat scattering using the distribution of hollow spheres (DHS) method, to generate transit light curves between 0.45 and 2.5 $\mu$m. Results. We show that the transit depth is wavelength independent for optically thick tails, potentially explaining why only some observations indicate a wavelength dependence. From the 3D nature of our simulated tails, we show that their transit cross-sections are related to the component of particle ejection velocity perpendicular to the planet's orbital plane and use this to derive a minimum ejection velocity of 1.2 kms$^{-1}$. To fit the average transit depth of Kepler-1520 b of 0.87%, we require a high dust mas-loss rate of 7 $-$ 80 M$_\oplus$ Gyr$^{-1}$ which implies planet lifetimes that may be inconsistent with the observed sample. Therefore, these mass-loss rates should be considered to be upper limits.Comment: 22 pages, 22 figures, accepted for publication in A&

Resolving the compact dusty discs around binary post-AGB stars using N-band interferometry

We present the first mid-IR long baseline interferometric observations of the circumstellar matter around binary post-AGB stars. Two objects, SX Cen and HD 52961, were observed using the VLTI/MIDI instrument during Science Demonstration Time. Both objects are known binaries for which a stable circumbinary disc is proposed to explain the SED characteristics. This is corroborated by our N-band spectrum showing a crystallinity fraction of more than 50 % for both objects, pointing to a stable environment where dust processing can occur. Surprisingly, the dust surrounding SX Cen is not resolved in the interferometric observations providing an upper limit of 11 mas (or 18 AU at the distance of this object) on the diameter of the dust emission. This confirms the very compact nature of its circumstellar environment. The dust emission around HD 52961 originates from a very small but resolved region, estimated to be ~ 35 mas at 8 micron and ~ 55 mas at 13 micron. These results confirm the disc interpretation of the SED of both stars. In HD 52961, the dust is not homogeneous in its chemical composition: the crystallinity is clearly concentrated in the hotter inner region. Whether this is a result of the formation process of the disc, or due to annealing during the long storage time in the disc is not clear.Comment: 12 pages, 10 figures, accepted for publication in A &

SPITZER survey of dust grain processing in stable discs around binary post-AGB stars

Aims: We investigate the mineralogy and dust processing in the circumbinary discs of binary post-AGB stars using high-resolution TIMMI2 and SPITZER infrared spectra. Methods: We perform a full spectral fitting to the infrared spectra using the most recent opacities of amorphous and crystalline dust species. This allows for the identification of the carriers of the different emission bands. Our fits also constrain the physical properties of different dust species and grain sizes responsible for the observed emission features. Results: In all stars the dust is oxygen-rich: amorphous and crystalline silicate dust species prevail and no features of a carbon-rich component can be found, the exception being EPLyr, where a mixed chemistry of both oxygen- and carbon-rich species is found. Our full spectral fitting indicates a high degree of dust grain processing. The mineralogy of our sample stars shows that the dust is constituted of irregularly shaped and relatively large grains, with typical grain sizes larger than 2 micron. The spectra of nearly all stars show a high degree of crystallinity, where magnesium-rich end members of olivine and pyroxene silicates dominate. Other dust features of e.g. silica or alumina are not present at detectable levels. Temperature estimates from our fitting routine show that a significant fraction of grains must be cool, significantly cooler than the glass temperature. This shows that radial mixing is very efficient is these discs and/or indicates different thermal conditions at grain formation. Our results show that strong grain processing is not limited to young stellar objects and that the physical processes occurring in the discs are very similar to those in protoplanetary discs.Comment: 22pages, 50 figures (in appendix), accepted for A&